Graphene nanoribbons on vicinal SiC surfaces by molecular beam epitaxy 
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We present a new method of producing a densely ordered array of epitaxial graphene nanoribbons 
(GNRs) using vicinal SiC surfaces as a template, which consist of ordered pairs of (0001) terraces 
and nanofacets. Controlled selective growth of graphene on approximately 10 nm wide of (0001) 
terraces with 10 nm spatial intervals allows GNR formation. By selecting the vicinal direction of 
SiC substrate, [1100], well-ordered GNRs with predominantly armchair edges are obtained. These 
structures, the high density GNRs, enable us to observe the electronic structure at K-points by 
angle-resolved photoemission spectroscopy, showing clear band-gap opening of at least 0.14 eV. 

PACS numbers: 73.22.Pr, 68.37.-d, 68.65.-k, 79.60.Jv 
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Graphene nanoribbons (GNRs) are attracting increas- 
ing attention in nanoelectronic applications and solid 
state physics, where the band-gap opening or modifica- 
tion of the electronic structure at K-points is a central 
interest^. The electronic structure at K-points in GNRs 
has theoretically^"'^ and experimentally^"^ ^ been demon- 
strated to depend on the type of edge geometry: armchair 
or zigzag. Semiconducting characteristics arc expected in 
the case of armchair edges owing to the band-gap open- 
ing at K-points^. As the width of GNRs is reduced, the 
gap is increased by both electron confinement and edge 
effects. However, realization of GNRs with atomically 
well-defined edges and providing experimental evidences 
of the gap opening at K-points remain challenging. Here 
we demonstrate a new approach for producing a densely 
ordered array of ahgned GNRs, which is advantageous 
for investigating physical properties macroscopically, on 
unique SiC surfaces as templates via molecular beam epi- 
taxy (MBE) and show band-gap openings at K-points 
visualized by angle-resolved photoemission spectroscopy 
(ARPES). Fabrication of GNRs by a conventional litho- 
graphic technique was first reported by Chen^ and Han^ . 
GNRs as narrow as 15'--^20 nm were obtained and the 
band-gap openings were shown; the band-gap value was 
dependent on the ribbon width. Since then, new and 
more advanced approaches have been proposed for GNR 
fabrication. These include unzipping of carbon nan- 
otubes (CNTs)'', chemical®, sonochemical^, bottom-up 
molecular precursor methods^*^ and selective thermal de- 
composition at facets on SiC substrates^^. Each method 
is unique and all have advantages and disadvantages. The 
unzipping of CNTs and chemical methods are appropri- 
ate for mass production but in their present state, re- 
quire highly complicated treatment and lack sufficient 
quality control. The primary advantage of the molecu- 
lar precursor method is that it enables very precise con- 
trol of the width and edge structure; however, so far, 
it only works on metal surfaces. The thermal decom- 



position of SiC substrates, utilized in reference 11, is 
now recognized powerful method to obtain high- 

quality graphene with few structural defects^^ owing to 
the epitaxial characteristic of the initial carbon layer 
with (6a/3 x6v^)R30°(6R3) registry on SiC(lx l)".i4. 
GNRs, approximately 40 nm wide, were achieved by the 
thermal decomposition method with selective growth on 
facets created by reactive ion-etching on a SiC substrate. 
However, in this method, C atoms are supplied inter- 
nally from SiC surface steps, and accordingly, unstable 
growth^^"^^ and morphological ffuctuation^^ are induced. 
Although all these methods have succeeded in forming 
GNRs of different quality, the physical characteristics of 
the resulting GNRs, in particular evident and direct ob- 
servations of band-gap opening, have not been demon- 
strated. 

Here we report a new method for the fabrication of 
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FIG. 1. An AFM image and a structural model of SiC 
nanosurface. (a) An AFM 3D view of the vicinal SiC surface 
after hydrogen etching, showing periodic array of terraces and 
facet structures, (b) A schematic drawing of the SiC surface 
in (a). Each pair of array consists of the (0001) terrace and 
(llOn) facet and is ordered with the periodic distance of ~20 
nm. The width of each (0001) terrace is ~10 nm. 
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GNRs on a SiC substrate by molecular beam epitaxy 
(MBE). C atoms are supplied externally, which effec- 
tively addresses the problems described in reference 11. 
Growth selectivity is achieved through the preferential 
nucleation of GNRs on terraces of vicinal SiC surfaces. 
The off-axis SiC substrate, often used in the fabrication of 
SiC electronic devices, is especially noted here because it 
shows a unique periodic surface structure. We have stud- 
ied off-axis (vicinal) SiC surfaces that are intentionally 
off-cut from a (0001) plane and found self-ordered peri- 
odic structures consisting of pairs of a (0001) basal plane 
terrace and a (iTOn) nanofacet (n = 35^37) with a char- 
acteristic periodicity of ~20 nm, as shown in Fig. 1. The 
structural model, illustrated in Fig. 1(b), was derived 
from high-resolution transmission electron microscopy 
images. The surface structure is formed by phase sep- 
aration, quantized step-bunching and ordering^®. Such a 
unique periodic SiC surface, hereafter SiC nanosurface, 
should be a good template for growing a massive array 
of GNRs. 

In the present study, vicinal 6H-SiC substrates (Si- 
face, 4° off toward [1100]) first underwent H2 gas etch- 
ing to produce a SiC nanosurface. The resultant SiC 
samples were then loaded into an ultra-high vacuum 
chamber, heated by direct current and cleaned under 
Si-flux at 1050°C. In situ reflection high-energy elec- 
tron diffraction (RHEED) indicated that the initial broad 
(a/3 X ■\/3)R30° (R3) structure of native oxide was changed 
to clear (3x3) of clean SiC(lxl) with Si adlayers^^, 
which are free from carbon and oxide contaminations. 
After terminating the Si-beam at the same temperature. 
Si adatoms R3 structure appeared owing to the evapo- 
ration of Si atoms in the Si adlayers^°. The Si adatom 
R3 surface is important to achieve high-quality GNRs on 
SiC by MBE, because previous highly Si-rich (3x3) struc- 
ture induces undesired SiC growth owing to residual Si 
atoms along with C atoms that are externally supplied, 
resulting in a collapse of the periodic morphology. Af- 
ter that, C atoms were supplied at a constant deposition 
rate by heating resistive carbon plates at 2200°C. The 
RHEED pattern was monitored in situ at 15 keV and 
used to terminate the growth. Three samples were pre- 
pared by changing the growth time: 5, 20 and 33 min. 
Each sample was evaluated by atomic force microscopy 
(AFM), RHEED, low-energy electron diffraction (LEED) 
and Raman spectroscopy. Finally, the samples were ex- 
posed to hydrogen at 600° C for 1 h to transform the 
6R3 structure into quasi-free-standing graphene by hy- 
drogen intercalation^^. In addition, scanning tunneling 
microscopy (STM), polarized Raman spectroscopy and 
ARPES were performed. 

Figure 2(a) shows the RHEED image after 20 min 
growth. The electron beam was irradiated parallel to the 
vicinal direction, i.e. perpendicular to step edges [1100]. 
Three sets of diffraction patterns are visible in this im- 
age: xl SiC surface, x3 due to R3, and satellites due 
to 6R3. Figure 2(b) shows the evolution of RHEED in- 
tensity with growth time at selected positions; namely. 
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FIG. 2. In situ RHEED analysis during MBE growth and 
a LEED image after growth, (a) A RHEED image after 20 
min growth, where both R3 and 6R3 structures are visible, 
(b) A LEED image after 33 min growth, indicating clear 6R3 
satellite spots, (c) RHEED intensity profiles as a function of 
growth time. The intensities of a specular spot, an R3 and a 
6R3 streak, as indicated in (a), are monitored. Note that the 
initial R3 intensity is increased after 5 min growth and the 
R3 and 6R3 intensities show a complementary relationship. 



a specular spot, R3 and 6R3 streaks, as indicated in 
Fig. 2(a). The diffraction intensity of the R3 structure 
(2/3 streak) increased initially and decreased after ~5 
min. The initial increase in the intensity is probably re- 
lated to the structural change to the C-ricli R3 surface. 
The detailed mechanism of structural transformation at 
this stage is not yet clear. We can speculate that the 
C atoms may incorporate into the Si-rich phase^^ while 
maintaining the same surface geometry of R3 but result- 
ing in higher intensity of the R3 streaks. This struc- 
ture is now under investigation by in situ LEED I-V 
analysis. As the R3 intensity started to decrease, 6R3 
streaks appeared and increased in intensity. This com- 
plementary relationship between the intensity proflle of 
R3 and 6R3 structures implies structural transformation 
from R3 to 6R3 by the incorporation of C atoms. Af- 
ter 33 min growth, the transformation was complete; no 
R3 signal was visible. At this stage, the 6R3 layer cov- 
ered the (0001) terraces, as shown in the LEED image in 
Fig. 2(b). Note that no additional (overlap) spots except 
for SiC and the 6R3 structure are seen in this image, in- 
dicating that growth takes place selectively on (0001) ter- 
races. The Raman spectrum of the sample grown for 33 
min typically indicates no G' (2D) signal [see Fig. 4(a)], 
which also supports the formation of the 6R3 structure. 

Next, each sample was hydrogen intercalated to ob- 
tain quasi- free-standing graphene on the SiC surface (21). 
LEED (not shown) clearly indicated the absence of 6R3 
satellites due to structural transformation and shows only 
superposition of clear (1x1) SiC and (1x1) graphene 
patterns. In the AFM height-mode, a slight contrast 
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FIG. 3. AFM analyses of GNRs. AFM height (upper) and 
phase images (lower) of the samples after 20 min growth (a 
and c) and 33 min growth (b and d) and the hydrogen in- 
tercalation, (e) represents the cross-section at the line A-B 
indicated in (a) and (c). ~0.25 nm higher regions shaded in 
(e) show darker contrast in the phase image, i.e. graphene on 
the hydrogen-terminated SiC surface. 



was visible on the (0001) terraces of the sample after 
20 min growth, as shown in Fig. 3(a). No contrast was 
observed before the hydrogen treatment. This contrast 
results from the difference in height between graphene 
(previously 6R3) and the R3 regions. As the 6R3 area 
transformed to graphene, its height increased owing to 
hydrogen intercalation at the interface, whereas the R3 
area was probably etched in the hydrogenation process. 
The height difference is ~0.25 nm, which is in good agree- 
ment with the cross-sectional height of graphene on the 
H-tcrminated SiC surface^"^. This is seen more clearly in 
Fig. 3(c), which shows the AFM phase image contrast for 
the corresponding area^'' and is also evident from a com- 
parison of cross-sectional height and phase profiles at the 
line A-B in Figs. 3(a) and (c). The sample grown for 33 
min, i.e. covered with GNRs, shows a surface morphol- 
ogy that is very similar to that of the SiC nanosurface, as 
shown in Fig. 3(b), but areas of dark contrast at terraces 
are evident in the phase image, as shown in Fig. 3(d), 
owing to the formation of graphene. 

The Raman spectra of the GNR sample (33 min 
growth) before and after hydrogen intercalation are 
shown in Fig. 4(a). The spectrum before hydrogen in- 
tercalation (as-grown) indicates typical features of the 
6R3 structure, including lack of D- and G'-bands be- 
cause of the absence of Dirac cones. In contrast, the 
GNR sample (after hydrogen intercalation) reveals evi- 
dent D-, D'-, G-, and G'-bands. The D- and G'-bands 
show a fairly narrow peak fitted with a single Lorentzian 
function. The full-widths at half-maximum of the D and 
G' peaks are 31 and 49 cm~^, respectively. The G-band 
peak at 1583 cm^^ is due to LO phonons at armchair 
edges, which are softened owing to the Kohn anomaly 
effect^^. The deconvolution analysis of the G-band in- 
dicates that there is a small additional peak at a wave- 
length 15 cm~^ higher, which is probably due to either 
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FIG. 4. A Raman analyses of GNRs. (a) Raman spectra of 
the sample after 33 min growth before (blue) and after (red) 
hydrogen intercalation, (b) Polarization angle dependence of 
the D-band intensity. D-band intensity at each polarization 
angle 0, defined in the insert, is fitted as a function of cos*S 
shown by a dotted curve. 



bulk LO phonons or TO phonons at armchair edges in 
the GNRs. The relatively strong D-band intensity due to 
the presence of high-density armchair edges and possible 
point defects should be noted. The edge character can 
be confirmed by the polarization angle dependence of the 
incident laser on the D-band intensity^^. The integrated 
intensity of the D-band at polarization angle 6 is plotted 
in Fig. 4(b). These are fitted as a function of cos^6' , 
suggesting predominantly armchair-type edges. This is 
also supported by the STM observation; line nodes^^ in 
the vicinity of the edges due to the interference effects of 
wave functions were recognized. 

The band structures of the GNR samples were stud- 
ied using ARPES. In Fig. 5(a), the ARPES result along 
the F-K-M lines of the graphene surface Brillouin zone 
is shown as a photoemission intensity image. A linear 
valence band dispersion verifies the presence of a single 
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FIG. 5. Intensity map and EDCs of the ARPES spectra 
around the K-point. Intensity map (a) and EDCs (b) of the 
ARPES spectra around the K-point of the samples after 33 
min growth. The spectra are taken along the F-K-M line. 
The red dots in (a) indicate the positions of the EDO peaks. 
The intensity map clearly indicates a linear dispersion of the 
valence band of the single-layer graphene. The EDCs show 
the folding of the valence band at the K-point below Ef. 
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layer graphenc. Note that a conduction band is invisible 
and the band is folded at the K-point shown in 

Fig. 5(b). The red dots in Fig. 5(a) indicate the peaks of 
the energy distribution curves (EDC) shown in Fig. 5(b) 
and represent band dispersion around the K-point. No 
states are detected between the Fermi energy (Ep) and 
valence band maximum, and a band-gap is opened at 
the K-point. The minimum band-gap of our GNRs can 
be 0.14 eV, assuming the conduction band minimum at 
Fip. However, based on the results obtained in the hy- 
drogen intercalation of the 6R3 sheet, the band-gap value 
is expected to be more than twice of 0.14 eV. A quasi- 
free-standing monolayer graphene shows a slight p-type 
doping nature^^'^*, i.e. the Dirac point should be above 

In summary, by templating vicinal SiC surfaces con- 
sisting of ordered pairs of (0001) terraces and nanofaccts, 
GNRs are formed by MBE. The carbon atoms supplied to 
such a surface selectively organize a graphene network on 
(0001) terraces and form an epitaxial 6R3 layer. Hydro- 
gen intercalation results in transformation to quasi-free- 
standing graphene and thus GNRs. The edge character- 
istic of the GNRs grown under these experimental condi- 



tions indicates armchair-type edges. The massive arrays 
of GNRs, ^10 nm in width, indicate apparent band-gap 
openings of at least 0.14 eV at the K-point observed by 
ARPES. A larger band-gap energy should be expected 
in GNRs ~^5 nm in width. Such GNRs can be grown on 
vicinal 4H-SiC surfaces whose terrace width is ~5 nm^^. 
These structures arc also interesting in terms of catalytic 
phenomena, which require high density of edges to de- 
tect. 
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